Spin Relaxation in a Quantum Dot due to Nyquist Noise 



O 

o 

(N 

o 

Q 



Florian MarquardtQ and Veniamin A. Abalmassov^j] 

1 Department of Physics, Yale University, New Haven CT, 06520, USA and 
2 Institute of Semiconductor Physics SB RAS and Novosibirsk State University, 630090 Novosibirsk, Russia 

(Dated: 2nd February 2008) 

We calculate electron and nuclear spin relaxation rates in a quantum dot due to the combined 
action of Nyquist noise and electron-nuclei hyperfine or spin-orbit interactions. The relaxation rate 
is linear in the resistance of the gate circuit and, in the case of spin-orbit interaction, it depends 
essentially on the orientations of both the static magnetic field and the fluctuating electric field, as 
well as on the ratio between Rashba and Dresselhaus interaction constants. We provide numerical 
estimates of the relaxation rate for typical system parameters, compare our results with other, 
previously discussed mechanisms, and show that the Nyquist mechanism can have an appreciable 
effect for experimentally relevant systems. 
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I. INTRODUCTION 

The dynamics of electron and nuclear spins in quantum 
dots has become a focus of research recently, with part of 
the motivation stemming from the prospect of potential 
applications in the context of spintronicsi or quantum 
computation^. Of particular importance are relaxation 
processes leading to decoherence, which have received 
much theoretical 4 ! 5 ! 6 : 7 : 8 : 9 and experimental 10 ! 11 ! 12 : 13 : 14 : 15 
attention within the last few years. 

Relaxation between Zeeman-split spin levels involves 
energy dissipation, which, in quantum dots with their 
discrete orbital spectrum, requires coupling to external 
degrees of freedom (a bath). In most theoretical works, 
phonons have played the role of such a bath 4 * 5 * 6 .. Here, 
we will consider the Nyquist fluctuations of electric fields 
produced by nearby gate electrodes as an alternative 
source of dissipation. Decoherence due to Nyquist noise 
has so far been considered extensively only for charge- 
based quantum computation proposal o 16 ! 17 ' 18 . The po- 
tential contribution to the spin relaxation rate of mag- 
netic field fluctuations (due to Nyquist current noise) has 
been found to be very small^. More recently, decoher- 
ence of phosphorous electron and nuclear spins in silicon 
due to Nyquist-noise induced fluctuations of the hyper- 
fine constant has been analyzed in Refs. 00- 

In the present work, we will calculate the electron spin 
relaxation rate in a single lateral quantum dot (Fig. |T|> . 
due to Nyquist fluctuations of a gate voltage, combined 
either with electron-nuclei hyperfine interaction or spin- 
orbit interactions. In the case of hyperfine coupling, this 
also contributes to nuclear spin relaxation within the 
quantum dot. In the case of spin-orbit coupling, the spin 
relaxation rate displays a striking dependence on the di- 
rections of both the fluctuating electric field as well as 
the static magnetic field. The spin-orbit Nyquist mech- 
anism considered in this work can become as efficient as 
coupling to piezoelectric phonons (the most important 
other mechanism) in realistic experimental setups. The 
rate grows linearly in resistance of the gate circuit. 

In the following, we will first introduce our model 




Figure 1: (Color online) Typical gated lateral quantum dot 
structure considered in the text, with gates connected to an 
electrical circuit of impedance matrix Z. The fluctuating elec- 
tric field produced by one particular gate is shown schemati- 
cally. 
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Figure 2: (a) Equivalent circuit scheme used in our model; (b) 
Energy diagram, showing the second-order relaxation process 
releasing the Zeeman energy hu) into the electromagnetic en- 
vironment. 



Hamiltonian (Sec. [H]), then describe the general cal- 
culation of spin flip rates in second order Golden Rule 
(Sec. IIII|I . and apply this scheme to both the hyper- 
fine (Sec. IIV|I and spin-orbit (Sec. |Vj) interactions. We 
will comment on the relation to recent spin relaxation 
measurements 1 ^ 1 ^. Finally, we will explain how the re- 
sults may be applied to arbitrary gate geometries and 
impedances (Sec. IVIfl and comment on deviations due to 
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other dot potentials or spatially inhomogeneous electric 
field fluctuations (Sec. IV 11(1 . 



II. THE MODEL 

We start by considering a single electron in a quan- 
tum dot formed in a 2DEG by a circularly symmetric 
parabolic lateral confining potential, subject to a homo- 
geneous static magnetic field B: 



Ha 



TV 

2m 



mQ 2 , 



p 2 + V(z)+gp B S-B, 



(1) 



where iv = p (q/c)A is the kinetic electron momentum, 
q = — |e| the electron charge, p the in-plane position 
vector, m and g the effective electron mass and g-factor 
respectively, ps > the Bohr magneton, V(z) the trans- 
verse confining potential, and f2o the lateral frequency. 

The solutions of the in-plane orbital part are Fock- 
Darwin states with energies n + fno + + 

where lo± = ^/ S1q + (uj c /2) 2 ± lu c /2 = ft ± uj c /2 and the 
cyclotron frequency is uj c = \q\B z /(mc). We will need 
the orbital ground state <^oo and the wave functions of 
the first excited doublet, 4> 00 = (7r/) _1 / 2 exp [— p 2 /(2l 2 )] 
and </>io = </>oi = (pe ilp /l)4>oo, with x + iy = pe tv and 
I = ^h/{mVL). 

We now consider spin relaxation due to the Nyquist 
fluctuations of an electric field acting onto the electron, 
being produced by a nearby gate (see Figs.n an d|21). We 
will at first assume the direction of the fluctuating field 
to be fixed (Fig. EJ) , and later describe how the analysis 
may be applied to arbitrary realistic gate geometries (Sec. 
IVI|) . This interaction is described by 



V E = qEp E = qE(x cos ( + y sin Q , 



(2) 



where E is the electric field amplitude, and pe denotes 
the projection of the electron's coordinate along the di- 
rection C of the in-plane field. We have used a dipole- 
approximation, assuming the electric field to be approxi- 
mately constant across the dot. This purely orbital inter- 
action can lead to spin decay only when combined with 
some spin-dependent part t4 p i n of the Hamiltonian (to 
be specified below). Thus, the total perturbation added 
to Hq reads: 



SH = V n 



V E . 



(3) 



Besides, the full Hamiltonian also contains a term de- 
scribing the electromagnetic environment, which deter- 
mines the dynamics of the electric field. 



III. SPIN FLIP RATE IN SECOND ORDER 
GOLDEN RULE 

According to Fermi's Golden Rule, the perturbation 
Ve produces transitions at the rate 

T = (2rr/n 2 )\(f'\p E \t')\ 2 e 2 (EE) u . (4) 



Here Tuo = Ei> — Efi is the energy absorbed by the en- 
vironment, and (EE) U = (2-k)- 1 J dt e iuJt (E(t)E) is the 
spectrum of the electric field. The spectrum may be re- 
lated to the impedance of the gate circuit (Fig. [5J by 
using the Fluctuation-Dissipation Theorem (FDT)i£: 



(5) 



Here d is an effective distance determining the conversion 
between voltage drop and electric field (distance between 
capacitor plates in the simplest case; see Sec. IVII for the 
general case). For the circuit of Fig. |3 the effective 
resistance responsible for the voltage fluctuations is given 
by R{uj) = R/(l + {RCu) 2 ). 

Spin-flip transitions are possible only because the 
matrix-element in Eq. Q is to be evaluated in pres- 
ence of the spin-dependent part V^pin, leading to a finite 
admixture of different spin states. According to standard 
stationary perturbation theory (leading to second order 
Golden Rule), we have: 

(/'I Pe |0 - (/II p E G(£ n )V spin + V spiB G(£ fl )p E |iT>. 

Here, the resolvent for H is G{£) = {£ — H ) For the 
transition between the Zeeman sublevels of the orbital 
ground state, we have £^ = Sq + fku/2 = £f± + Hui, where 
huj = \gpsB\ is the Zeeman energy, and en is the energy 
of the orbital ground state. Eqs. Q and iJBJ describe a 
second-order transition from to |/|). Note that pe 
connects the lateral orbital ground state c/>oo to <pio and 
0oi only. 

We will now consider two specific spin-flip perturba- 
tions V SDin . 



IV. HYPERFINE INTERACTION 

The hyperfine contact interaction of an electron spin 
S with nuclear spins Ij at positions Rj has the form |20| 



7 H F = y>oA i S-I J - £(r-R,-), 



(7) 



where the Aj are the hyperfine coupling constants (de- 
pending on species), and vq is the volume of the unit cell. 
We will neglect the nuclear Zeeman splitting. 

In the limit of small Zeeman splitting, uj ^ Qq, we 
obtain for the amplitude (JSJ): 



(/'I Pe |0 



HF 



-2lv - 



-^ 2 /(2»o) 2 ) 1/2 



(8) 



The initial and final states of the nuclear spin 
system are denoted in, Jn- We have defined 
H z j>Pj) = xl(zj)^oo(PjMw(Pj)\, for the nuclei posi- 
tions (pj, (fj, Zj), and Xo(z) is the ground state of motion 
along the z-direction. 
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Inserting into Eq. {2J and averaging over an ensemble 
of unpolarized uncorrelated nuclear spins, we obtain the 
following expression for the electron spin relaxation rate 
due to Nyquist noise and hyperfine coupling: 

r _4tt A 2 v r, R(u)u(l+u; 2 /(2n ) 2 ) 

HF 3 l{i + I) ^ m ^ 2 Z()d2 Rq 1 _ e _ n u/ {kB T) > 

(9) 

where A 2 — ^2jA 2 with summation over all nuclei in 
the unit cell, Rq — h/e 2 is the quantum of resistance, 
•q — zq1 2 J^>(z, p) 2 dzdp is a dimensionless form factor in- 
volving an average over nuclei positions (z,p,ip), and zq 
is the 2DEG thickness. 

In the following, we provide numerical estimates using 
typical parameters for GaAs QDs: hyperfine coupling 
constant 21 A 2 ~ 1.2 x 10 -3 meV 2 , nuclear spin value / = 
3/2, unit cell volume v = (5.65 A) 3 , effective electron 
mass m = 0.067m , \gfi B \ = 0.025 meV T" 1 , hu c /B z = 
1.76 meV T _1 , zq = 10 nm, and a geometrical factor 
of rj = 9/(16 • 47r) for the approximate solution of an 
inversion layer (triangular well) potential (see, e.g., Ref. 

Thus, the numerical value of the electron spin- flip rate 



HF 



0.6Hz x «(W(lfE£Vy B_ 
Rq \ d J \ hfl J IT 



1 +0.8 x 



/ 1 mcV B z \ ' 

\ nn it J 



1 -exp [-0.3 x 



(10) 



Concerning R(cu) = R/(l+ (RCu) 2 ), a lower bound for 
the cutoff frequency 1 / (RC) is determined by the charg- 
ing energy E c of the dot (involving the total capacitance 
Ctot > C), in the form (E c /H)(Rq/R). For current 
GaAs experiments, a typical value of E c = 1 meV yields 
RCu ~ 10- 1 (R/R Q )(B/1T), such that R(oj) deviates 
from R only at rather high magnetic fields, as long as 
R < Rq. As a function of R, the maximum relaxation 
rate is reached at R = \/{Cuj). In typical GaAs quantum 
dots this corresponds to R max = 4QRq(1T/B), 

For a confinement frequency of HCIq — 1 meV, the rate 
is equal to about 0.1 Hz at T = 1 K and B = 0.1 T, 
with R/Rq = 10~ 2 , d = 0.5 /im. This is comparable 
to electron spin relaxation due to hyperfine interaction 
involving piezoelectric phonons^ at the same values of 
Oo, T and B. Since the Nyquist mechanism considered 
here yields a zero-temperature rate linear in transition 
frequency (instead of cubic as in Ref. ®, it dominates at 
smaller magnetic fields B < 0.1 T. The relaxation rate 
can reach large values if the resistance approaches i? max ; 
e.g. at T = 1 K and B = 1 T (with R max = 40R Q ), we 
would have Tjjf — 0.3 kHz, much larger than the rate 
due to electron-phonon coupling 5 . 

If new electrons are supplied to the dot in a trans- 
port situation or an additional electron spin relaxation 
mechanism is effective, then the present mechanism may 
also result in nuclear spin relaxation, with a rate equal 



to Thf/A^, with TV the number of nuclei. For the param- 
eters used here (with N ~ 5 x 10 5 ), and at T = 1 K, 
R = 10~ 2 Rq, the nuclear relaxation rate turns out to 
be negligibly small (10 -7 Hz). In comparison, a recent 
transport experiment in GaAs quantum dotsi 3 - found nu- 
clear spin relaxation times of the order of 10 minutes at 
T = 100 mK and B — 40 mT (cf. Ref. for relevant 
theory). In order for the present Nyquist mechanism to 
yield a comparable rate, a resistance R ~ i? ma x ~ 10 3 i?g 
would be required. 



V. SPIN-ORBIT INTERACTION 

The other, numerically much more important spin- 
relaxation mechanism we will consider is due to the com- 
bination of Nyquist noise and spin-orbit coupling: 

VSO = Ol(TT x &y ~ TTy(T x ) + P{TT y d y - 7t x & x ). (11) 

The first term in 111|) is the Rashba term arising from 
any structural inversion asymmetry of the heterostruc- 
ture, the second is the (linear) Dresselhaus term, a con- 
sequence of the bulk inversion asymmetry of the semi- 
conductor material. The main crystallographic axes are 
assumed to be aligned with x, y, z. 

By means of the commutation relations (iK/m)it = 
[p, Hq] and [Hq, <x] = ig/is<r x B, we get: 



(12) 



where we have introduced the vector £ = (m/h)(—Px — 
ay, ax + (3y, 0). 

The contribution of the first term in Eq. (|12|) to the 
amplitude, Eq. ©, yields [/je,£<t] = 0. Thus, we obtain, 
in the limit lu <^ f2o : 



(f'\pE\i%o = gMB(/ll*[B x n]|iD, (13) 

Here, we have defined the vector ft = (cos(£ — 
7), — sin(C + 7), 0), where 7 depends on the ratio of spin- 
orbit coupling parameters: tan7 = a/ft. 

We parameterize the direction of the magnetic field as 
B = B(smd cosLp,sm9smip,cosd). Then the electron 
spin relaxation rate due to spin-orbit coupling combined 
with Nyquist noise is given by: 

a 2 +p 2 R{u) ^ J(c,7;0,y) hA] 



d 2 



R Q 0% 1 - e -^/( fe « T ) ' 



where 



I(C,rA<P) = (sin 2 (C + 7 )+cos 2 (C-7)) 

x (cos 2 6 + sin 2 6 sin 2 (ip - ip*)) (15) 

yields the angular dependence of the relaxation rate 
(see Fig. |3J), with the angle <p* defined by the condition 
tan (p* = — sin(£ + 7) / cos(( — 7) . 
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Figure 3: (Color online) Angular dependence of the spin re- 
laxation rate Fso on magnetic field direction (9, ip), according 
to Eq. 115II . for tp* — — tt/4. The shape will be rotated around 
the z-axis for different values of tp* , which is determined by 
the electric field direction £ and the spin-orbit coupling con- 
stant ratio expressed via 7. 



We note that, for fixed direction of the fluctuating elec- 
tric field, / vanishes for in-plane magnetic fields pointing 
into the direction tp* , i.e. /(£, 7 ; tt/2, cp*) = 0, V£,7 
(when £1 \\ B). In the particular case of equal spin-orbit 
coupling constants, a = ±/3 and thus 7 = ±7r/4, the first 
factor in / reduces to 2 sin 2 (£ ± 7r/4). This means the re- 
laxation rate may even vanish regardless of magnetic field 
direction, provided the electric field across the dot points 
along C = TV 4 : thus /(tt/4, ±tt/4 ; 9, tp) = 0, V6,(p 
(since Q = in Eq. Q13[0. In the case of an arbitrary 
gate geometry (treated in Sec. IVI|) . the rate usually does 
not vanish any more but can still have some pronounced 
directional dependence. 

Conversely, for a = ±/3 the rate may also vanish 
for arbitrary orientation £ of the electric field, if the 
magnetic field lies in the plane and tp = T-tt/4 , i.e. 
BC, ±7r/4; 7r/2, =ftt/4) = 0, VC- As pointed out in Ref. 
l24l for a = ±/3 the electron spin component a x =p a y 
commutes with the spin-orbit interaction Vsoj an d for 
B || (l,=pl,0) it commutes with Ho as well. Thus, the 
suppression of spin relaxation in this case is exact (see 
also Ref. 0). 

The numerical value for the relaxation rate is: 



r so « 0.4 kHz x 



R{lo) f l/im\ 2 / lmeV\ 4 



r q \ d J \ nn 



2q 



B 
FT 



/ l/xm\ ' 



1 — exp 



B IK 

-0.3 x 

IT T 



(16) 



with 1/A| = (m/h) 2 {a 2 + (3 2 ). 

The rate T so is about 1 Hz at T = 1 K, for B = 0.1 T, 
using the same parameters as for hyperfine interaction 
(and Aso = 1/ini, /(£, j;8,tp) = 1). At these parameters, 
this is comparable with the relaxation rate due to the 
combined action of spin-orbit and piezoelectric electron- 
phonon interactions 4,6 , which, however, vanishes like B°, 



instead of the B 3 dependence we have found here, such 
that the Nyquist mechanism dominates for lower fields. 
While a B 3 -dependence was also found in Ref. for a 
two-phonon relaxation process, the efficiency of that pro- 
cess decreases drastically when the temperature becomes 
lower than about 1 K. 

Very recently, the spin relaxation time in a one-electron 
quantum dot has been measured directly, using a time- 
resolved single-spin detection setupii Previously, it had 
been possible only to obtain lower bounds on this time 
(see Ref. ^3 for an example). At an external mag- 
netic field of 8T, the relaxation time was found to be 
on the order of 1msec. Using parameters M7o = 1 meV, 
d = 0.5 fj,m, and Aso = 1 (un, we indeed can obtain a 
rate of about 10 Hz by assuming a reasonable value 
of R/Rq = 10~ 2 (lower values of R are needed if the 
effective d is smaller). The observed magnetic field de- 
pendence is not inconsistent with a B 3 contribution, al- 
though more data are needed in this regard. Thus, the 
present Nyquist mechanism may be about as efficient as 
the piezoelectric mechanism of Ref. (i. for experimen- 
tally relevant situations, and it would be interesting to 
see whether the two effects can be distinguished in fur- 
ther measurements on similar setups (see below). 

In contrast, we do not expect Nyquist noise to con- 
tribute to spin relaxation rates in self-assembled quan- 
tum dotsi^, which are apparently well explained by piezo- 
electric electron-phonon coupling^ alone. 

As the rates Thf and Tgo depend on the relevant cir- 
cuit resistance R, it might be interesting to have a device 
where R can be controlled (e.g. in a superconducting gate 
circuit, with a strong decrease in R(lu) for T < T c and 
tunability by a magnetic flux; or using varactor diodes as 
tunable capacitors in an LC-circuit) . Inserting controlled 
dephasing by a large R could be a way to cross-check 
measurements of spin relaxation rates (e.g. R = 4_Rq 
would yield an easily measurable rate Tgo 011 the order 
of 10 7 Hz, for the previous parameters). 



VI. ARBITRARY GATE GEOMETRIES 

We now explain how our results may be extended to ar- 
bitrary experimental configurations for gated lateral dot 
structures (see Fig. QJ. 

In general, the dot experiences a fluctuating electric 
field which is due to the fluctuating voltages Uj on the 
gates j = 1 . . . Nq ■ (We assume these voltages to be 
measured relative to ground) The relation between the 
in-plane electric field components (I = 1,2) at the center 
of the dot and the voltages is linear, 



Ei 



N G 



(17) 



The real-valued matrix M depends on the geometry of 
the gates and the dielectric substrate (it has to be ob- 
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Figure 4: Electron spin relaxation rate due to Nyquist noise, 
Eqs. (Eni and JTHJ, at Ml = 1 meV, R/Rq = 10~ 2 , d = 
0.5 nm, Aso = Ifim, ip = 0, 9 = tt/2, < = 0, 7 = -7r/4. At 
low fields, the hyperfine rate (HF) dominates over the spin 
orbit rate (SO). 



anisotropy will generally remain due to the anisotropy of 
the gate geometry (or the impedance matrix). 

However, even if the geometry were perfectly known, 
and the matrix M were calculated using a numerical 
solver of the Poisson equation, it is currently difficult to 
go beyond estimates for realistic setups (like, e. g., those 
used in the experiments by the Munich 1 — or Delfti2ii4 
groups). This is because the impedance matrix of the 
gate circuit (and thus the noise properties) have never 
been studied in detail at the high frequencies under con- 
sideration, and the nominal circuit diagrams are probably 
valid only at low frequencies. Nevertheless, our estimate 
for the effective resistance entering the Nyquist noise, 
R/Rq ~ 10~ 2 , seems to be roughly comparable to the 
actual numbers (Munich group 2 ^) or at least cannot be 
ruled out at present (Delft groupSi). In this context it 
should also be notedSS. that part of the circuit is at higher 
temperatures, possibly contributing more noise than es- 
timated using the base temperature of the setup. 



tained from an electrostatic simulation for specific ge- 
ometries). This matrix has the dimension of an inverse 
length, and roughly corresponds to the factor d -1 we in- 
troduced in Eq. (J5J for our simple model system (which 
can be described by Nq = 2 and M12 = —Mn in this 
scheme) . 

The gates are nodes in an electronic circuit that is 
described by an impedance matrix Zy. Thus, according 
to the FDT, we have: 

We now make the reasonable assumption that the mag- 
netic field is not strong enough to appreciably influence 
the impedance matrix of the circuit. Then the Onsager- 
Casimir relation Zjij(B,uj) = Zjji(—B,u>), evaluated 
at B = 0, implies reciprocity Zfj = Z^y and thus 
(Z rj (u) + Z* f (oj))/2 = ReZ rj (u). Inserting Eq. (TgJ 
into Eq. 1171) . we have 

= 1 ^/I/t E M Vj ,M lj K & Z j , j {u>) . (19) 
j'j 

Reciprocity makes the matrix (EiiEi) real- valued and 

symmetric. Physically, this means there is no circu- 
lar component to the thermal electric field fluctuations. 
Thus, at given lo, we can diagonalize this matrix by 
choosing two orthogonal directions in the 2DEG plane. 
The total spin relaxation rate becomes equal to the sum 
of the rates for these two directions, obtained according 
to our previous description. The angle C m our nota- 
tion denotes the direction of the given principal axis, and 
(^EE^ is equal to the corresponding eigenvalue. As a 

consequence, the total spin-orbit rate Tgo will n °t van- 
ish completely for any magnetic field direction, but an 



VII. OTHER DOT SHAPES AND BEYOND 
DIPOLE APPROXIMATION 

If we consider an anisotropic or, in general, non- 
parabolic lateral confinement potential, the orbital 
single-particle wave functions and energies will change. 
Most importantly, this results in nonvanishing transition 
matrix elements between the ground state and more than 
just two excited orbital states. The same is true if we go 
beyond the dipole approximation, i.e. we no longer as- 
sume the fluctuating electric field to be constant across 
the dot. 

Thus, the sum over intermediate states implicit in Eq. 
(|tj|) extends over more states. However, generally speak- 
ing the qualitative picture of our previous analysis does 
not change, and neither do the quantiative estimates, 
provided the deviations from the parabolic shape are not 
too large. This is because transitions via higher inter- 
mediate states are suppressed anyway by larger energy 
denominators, and, in particular, no extra channel via 
low-lying energies opens up. 



VIII. CONCLUSIONS 

In summary, we analyzed electron spin relaxation due 
to the combination of Nyquist noise with both hyper- 
fine and spin-orbit interaction. For the case of spin-orbit 
coupling, a remarkable dependence on the directions of 
magnetic and electric fields is observed. For equal spin- 
orbit coupling constants (realizable by tuning the Rashba 
term 24 ), the rate may vanish exactly for a particular in- 
plane magnetic field direction (independent of gate ge- 
ometry) or, for arbitrary angles of the magnetic field, 
become very small provided the fluctuating electric field 
points predominantly into a certain fixed direction. The 
contribution to the spin relaxation rate may be as large as 
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that of the most important other mechanism (piezoelec- 
tric electron-phonon coupling) , for relevant experimental 
parameters. Furthermore, if a larger (preferably tunable) 
gate resistance is relevant for a given experimental setup, 
or microwave noise is deliberately applied to the gate, 
the mechanism analyzed here may dominate over an ex- 
tended range of magnetic field values and might be more 
easily distinguished from other mechanisms. 
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